* label the additional variable ------------------------------------------------
label variable hdp_jun             "Pre-test"

* run the regressions models and save them to memory ---------------------------

eststo: reg hdp_nov treatment hdp_jun recruitment_pool, cluster(province)

eststo: reg hdp_nov treatment hdp_jun recruitment_pool                          ///
treatment_nonterror casualty_district kurdish_district akp_district             ///
higher_education min_margin turnout_jun, cluster(province)

eststo: reg hdp_nov treatment treatment_multi hdp_jun recruitment_pool,         ///
cluster(province)

eststo: reg hdp_nov treatment treatment_multi hdp_jun recruitment_pool          ///
treatment_nonterror casualty_district kurdish_district akp_district             ///
higher_education min_margin turnout_jun, cluster(province)

* print the results ------------------------------------------------------------

esttab,                                                                         ///
replace se b(3) se(3) stats(N vce r2) label nodepvars nonotes nomtitles         ///
order(treatment treatment_multi)                                                ///
substitute(vce "Clusters" cluster "81")                                         ///
title("Table S15: Regression models of HDP vote share")

* save the results -------------------------------------------------------------

esttab using "./tables/table_s15.tex",                                          ///
replace se b(3) se(3) stats(N vce r2) label nodepvars nonotes nomtitles         ///
order(treatment treatment_multi)                                                ///
substitute(vce "Clusters" cluster "81")                                         ///
title("Regression models of HDP vote share")

* remove the models from memory ------------------------------------------------

eststo clear